Abstract. We propose a three component model consist- 
ing of minihalos and galactic halos with embedded thin 
discs for absorbers producing all the observed classes of 
intervening quasar absorption line systems. We show that 
this model, based on CDM cosmology, can explain most 
of the observed statistical distributions of various types 
of absorption systems. Use of the Schechter luminosity 
function for absorbers, on the other hand, is consistent 
with the observations only if the number of galaxies was 
larger in the past and reduced with time due to merg- 
. ers. A strong chemical evolution in the halos of galaxies 
is indicated by the observed properties of CIV lines. We 
discuss our results in the light of the recent observations 
of the absorption line systems. 

^ ' Key words: quasars- absorption lines-galaxies- 
<— ^ . abundances, galaxies-halos — -disc 
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1. Introduction 

Quasar absorption lines provide unique, sensitive and di- 
rect probes of the gaseous content of galaxies and inter- 
galactic medium (IGM) over look back times comparable 
to a Hubble time. Studies of these lines promise to pro- 
vide insight into the processes that regulate the evolution 
of the IGM and the galaxies and the mechanisms that 
determine the physical conditions in discs and in clouds 
embedded in the halo. In a typical quasar spectra there 
are several populations of absorption lines and it has be- 
come apparent that these different populations are likely 
the absorption signatures of different regions and condi- 
tions within the galaxies and the IGM. 

Narrow, heavy element absorption systems (HEASs) in 
the quasar spectra are believed to be produced by gas in 
the halos of intervening galaxies. This was demonstrated 
by the high success rate of detection of emission at the 
rcdshift of the intermediate redshift (z ~ 0.5-0.8) Mgll ab- 
sorption systems (Bergeron 1988; Steidel 1993). High res- 
olution studies of these heavy element lines have revealed 
that the lines break up into several distinct components 
which are most likely produced by clouds in the halos of 
the absorbing galaxies (Srianand & Khare 1994). Most of 
the Lyman limit systems (LLSs), which have a neutral 
hydrogen column density > 2xl0 17 cm~ 2 are associated 
with the HEASs. 

Damped Lya systems (DLASs) are produced by con- 
densations of neutral gas having neutral hydrogen column 
density Nhi > 10 20 - 3 cm~ 2 and are easily detectable in 
the spectra of background quasars (Wolfe et al. 1986). 
These are associated with metal lines of low ionization 
and often high ionized species as well. Although the exact 
nature of the DLASs has yet to be understood it is gen- 
erally believed that they are produced by the progenitors 
of present day galaxies. This conclusion is based on the 
following facts (i) the mass contained in these systems at 
redshifts >2 seems to be comparable to the stellar mass at 
present and (ii) the kinematics of the gas in DLASs seems 
to indicate their origin in rotating discs. This view has 
been questioned recently by several authors in the light of 



some of the new data on these systems. We will discuss 
these later on in this paper. 

Lanzetta et al. (1995) and Chen et al. (1998), from a 
comparison of galaxy and absorber redshifts showed that 
a dominant fraction of low redshift (z<l) Lyman alpha 
absorption systems (LAASs) arise in extended halos of 
galaxies. This was also corroborated by the discovery of 
heavy elements (Cowie et al. 1995) at redshifts around 
3 and of clustering in these systems at redshifts between 
1.7 and 3.8 (Srianand & Khare 1994). Recently Chen et 
al. (2000) have considered the statistical properties of the 
high redshift galaxies in the Hubble deep field and the 
statistical properties of Ly a absorption systems. They 
conclude that at least 50% and possibly as high as 100% 
of observed LAASs arise in extended envelopes of galaxies 
at z>2. It thus appears that a substantial fraction of the 
LAASs may be produced by the outer parts of galactic 
halos rather than in the IGM as believed earlier (Sargent 
et al. 1980). 

During the past few years, there has been a maturing 
of attempts to understand the nature of these intervening 
population, evolving from rather simple idealized mod- 
els (e.g. pressure confined clouds, gravitationally confined 
"minihalos", freely expanding clouds etc) to sophisticated 
numerical simulations that are able to reproduce many 
of the features of the absorption spectra in the context 
of the evolution of structure (Chiba & Nath 1997; Das 
& Khare 1999; Charlton & Churchill 1996; Kauffmann 
1996; hereafter K96, Gncdin & Ostriker, 1996 etc.). Das 
& Khare (1999) have considered a two component (mini- 
halos and galactic halos) CDM model to understand the 
LAASs and associated CIV absorption lines. K96 consid- 
ered the formation of disc in CDM models and showed that 
the observed properties of the DLASs can be accounted 
for by these models. This pure disc model for DLASs was 
questioned by Prochaska & Wolfe (1997) on the basis of 
kinematic structure of the DLASs. Subsequently Jedamzik 
& Prochaska (1998) have shown that only a finely tuned 
disc model within the framework of CDM could be made 
marginally consistent with the damped Ly a observations. 

The purpose of this paper is to explore a consistent 
model of minihalos and of galaxies which are assumed to 
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have spherical and photo ionized gaseous clouds embedded 
in their halos with a disc at the centre, to better under- 
stand the nature of the LAASs, HEASs and DLASs in a 
unified manner. In section 2, we describe the basic model 
with the underlying hypothesis. We discuss the orienta- 
tion of the discs inside the halo in section 3. We explore 
the field of parameters and constrain these parameter us- 
ing observations in section 4. We present the results in 
section 5 and 6 before drawing our conclusions in section 
7. A discussion of the results in view of some of the recent 
observations is presented in section 8. 

2. The basic model 

The two component model was recently considered by Das 
& Khare (1999) to understand the observations of Lyman 
a forest lines and their accompanying C IV lines in the 
framework of the hierarchical structure formation model 
in CDM cosmology. This model consists of minihalos and 
galactic halos. These are dark matter halos having 5% of 
mass in the form of baryons. The Minihalos have circular 
velocities, V c between 15 and 55 km s _1 and consequently 
have smaller virial radius (see eqn 2.1). The lower limit of 
15 km s _1 is roughly the value below which the pertur- 
bations in the baryonic mass are suppressed while 55 km 
s^ 1 is roughly the value below which the photo ioniza- 
tion by the external UV background suppresses cooling 
and therefore the formation of clouds and stars. These are 
chemically enriched if at all, by an earlier generation of 
stars. 

Halos with circular velocities larger than 55 km s _1 arc 
called galactic halos. Upper limit on V c for galactic halos 
is taken to be 500 km s _1 as the results presented below 
do not change if the upper limit is increased above this 
value. The gas in these halos can cool due to higher gas 
densities in here (see below) and form clouds and stars. 
Thus these halos have clouds embedded in it which are 
in pressure equilibrium with the hot intercloud medium. 
Their number density is assumed to have the same ra- 
dial dependence as the total density. Here, we extend this 
model to include the presence of isothermal discs inside 
the galactic halos which form as the gas cools, collapses 
and forms stars (K96). It may be noted that all the halos 
in our model, being gravitationally confined are assumed 
to be spherical. There have been suggestions in the lit- 
erature of the filamentary nature of the Ly a absorbers 
(Charlton et al 1993). These were based upon the large 
transverse sizes of these absorbers (> 100 kpc) obtained 
from the occurrence of common absorption lines in the 
spectra of multiple images of the lensed QSOs (Smette 
1995). Similar sizes along the line of sight (los) (as would 
be required for the spherical absorbers) would require too 
low a value of the external (IGM) pressure. This is how- 
ever, not a problem in our model as the galactic halos 
(where the lines observed by Smette (1995) at low res- 
olution would have originated) have clouds with sizes of 



about a few kpc embedded in them. The common absorp- 
tion lines in the images of lensed QSOs may be produced 
by different clouds in the same halo. Thus the assumption 
of spherical halos does not violate the constraints from 
the lensed QSO data. We note that the lensing data in- 
dicate that if the los to one image intercepts an absorber 
the probability that the los to the other image intercept- 
ing the same absorber is > 90%. We will discuss this with 
reference to our model in section IV. 

Our model for the disc is based on the work of K96 ac- 
cording to which the gas collapses to a constant fraction 
of its initial virial radius to form a rotationally supported 
disc within the halos, the virial radius for spherical col- 
lapse model is given by 

r^O.lH-^l + zr^Vc (1) 
Here H Q is the Hubble constant. Using the results of N 
body simulation (Barnes & Efstathiou 1987) it can be 
shown that for a disc/halo mass ratio of 0.1 the equi- 
librium disc radius is ~ 10% of the virial radius of the 
halo. A spread around this factor is expected due to pos- 
sible spread in spin parameters of the dark halo (K96). We 
have, included an additional dependence of the disc size 
on redshift apart from the (l+z)~ 3 / 2 dependence which is 
due to the z dependence of r v ; r . This seems to be necessary 
to account for the observations as seen below. Thus the 
disc radius in our model is given by, 

R d = 0.1H- 1 (l+z) a - 3 / 2 V c , (2) 
where a is an unknown parameter whose value as well 
as the value of the thickness of the disc, t, is fixed by 
comparing the model results with observation of DLAs as 
explained below. 

Following K96, we assume all the discs to have flat ro- 
tation curves and a rotational velocity equal to the circular 
velocity of their surrounding dark matter halos. Assum- 
ing Kennicutt (1989) star formation law and a constant 
gas velocity dispersion of 6 km/s the gravitation stability 
condition gives the critical surface density as (K96) 



S crit (R) = O.59 V i k ^ 1) M 0P c- 2 



(3) 



R(kpc) 

where V is the velocity of the disc and R is the radial 
distance from the centre. The surface gas density in the 
disc remains close to the critical density with constant in 
fall, star formation and ejection of hot disc gas into the 
halo through supernovae explosions. So the surface gas 
density in the disc is also given by the above equation. 

The total hydrogen column density contributed by the 
disc for an impact parameter P for a face on disc can be 
written as, 

N ht(P) = (4) 

m being the average mass per gas particle which we as- 
sume to be mass of the proton, m p . The average number 
density, nn along a los through the disc is thus 

N£f(P) 



(5) 
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The total column density along a given los is the sum of 
the column density contributed by the disc, the column 
densities in all the clouds in the halo that lie along the los 
and the column density contributed by the hot medium 
along the los. We assume the number density of clouds at 
a radial distance r from the centre of the halo, rj c \(r) to be 

%i(0) 



*7d(r) 



1 + r 2 /r 



(6) 



r c being the core radius. The number of clouds per unit 
volume at the centre of the galaxy is given by 

^°) = f i7W (7) 

The value of f is fixed by comparing the model results with 
observations, as described below. 

The radial density distribution of the baryonic matter 
in the minihalos and in galactic halos is given by 



p(r 



(- 



-), 



(8) 



4 7r G (r 2 + r 2 ) v x v — arc tan x v 
where fi, is the fraction of baryonic matter in the halo, 
which we take to be 0.05 as mentioned above. r c is the 
core radius which is taken to be 10 and 100 kpc for mini 
and galactic halos respectively and x v = ^ ~ 15 (Chiba 
& Nath 1997). 



3. Orientation of the discs 

The orientation of the disc is an important factor while 
determining the column density contribution of the disc 
which is likely to increase due to the increased path length 
through an inclined disc. For the calculation of the contri- 
bution of the disc to the column density along a los with 
impact parameter P, we assume the disc to be randomly 
oriented with inclination i. At this orientation, the disc 
is projected on the plane of the sky on to an ellipse with 
area 7rRd 2 cos i. If the los is located at a random angle 9 
from the major axis of the ellipse then the radial position 
along the surface of the disc where it enters the plane of 
the disc is given by 



211/2 



(9) 



R in =P[cos 2 0+(^4) 

COS 1 

and the radial position along the surface of the disc where 
it leaves the plane of the disc is given by 

R out = [ P 2 cos 2 + ( P H£f + t tan i) 2 ] 1 / 2 . 

The column density for a particular ion along a los 
with impact parameter P, for minihalos is given by, 

Rh 

n H (r)fi „(r)rdr 



Ni on (P) 



Vr 2 - P 2 



(10) 



Here njj (r) is the number density of hydrogen, i?u is the 
radius of the halo and f; on (r) is the ionized fraction of the 
ion multiplied by Z, the abundance of the clement being 
considered. f; on depends on the ionization parameter T at 
the point under consideration. For galactic halos, when 



the los passes outside the disc the column density is given 

by 

Rh Rh 



,(P) = 2 J A(r)dr + 2 J B(r)dr. 



(11) 



p p 
Here, A(r) = ^H(r)W(r)gr)», cl (r)r dctcrmincs the contribu- 
tion from the galactic clouds, <r(r) being the cross section 
of the clouds. iVn(r) is the total hydrogen column density 
contributed by an individual cloud present at the radial 
distance r from the halo centre. B(r) = nH ^ fi °p^ r is the 
contribution to the column density from the intercloud 
medium. 

When the los with an impact parameter P passes 
through the disc, the column density for a particular ion 
is, 



N io „(P,M 



VR?„+(t/2) 2 

J A(r)dr + 



Rh 



Rh 

/ 



A(r)dr + 



y/Rf D +(t/2)2 



/ 

Rh 



B(r)dr 



V R out +(t/2) 2 



B ( r )dr + fi °" S(P) .(12) 



V R out +(t/ 2 ) 2 

The first two terms give the contribution of the clouds 
in the halo along the los, the next two terms give the 
contribution from the intercloud medium in the halo while 
the last term gives the contribution from the disc. The 
disc contribution to the column density along the sides of 
the projected ellipse is constant, but the halo contribution 
changes from point to point. Here [i = cosi and the factor 
of i accounts for the increased path length through the 

disc due to its inclination; fi °°^( p ) [ s the column density 
contributed by the disc in the face on position. 

The ionization parameter T at a radial distance r from 
the centre is given by T(r) = 4 * 10 , — ^ where an is 
the slope of the spectrum of the UV background radia- 
tion (assumed to be a power law). We assume o:q = -1.5. 
J _2i is the intensity of the UV background radiation in 
units of 10 -21 erg cm -2 s _1 Hz -1 str -1 at 1 Rydberg. c 
and h are the velocity of light and the Planck's constant 
respectively. The photo ionization code (cloudy 90) was 
run for a range of hydrogen column densities and a table 
^ vs. T was prepared for each value of hydrogen column 
density. These values are quite insensitive to the values 
of particle densities and abundances. This table was used 
to determine the value of for any value of T and Nh 
during the calculations. 

The column density distribution function for a partic- 
ular ion produced by halos with circular velocity Vi < 
V c < V u is given as, 
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d 2 N(z,N ior 
dN iorl dz 



H 



(1 + z) 



l,2 h 



(V c ,z)edV c 



dA(N ion , n) 
9N ion 1 



where |^^| 



v u 

= £(l + z) 1/2 | n(V c ,z)edV c 



1 A 7T 

: ° A(Nio.+AN^,M)-A(N iop , M ) and &rea 



A(Ni on , jLt) = J d# ^P 2 (8, n) is the projected area, for a 

0=0 

given value of i, inside which the column density is > N; on 
and n(V c ,z)dV c is the number of halos per unit volume 
with circular velocity between V c and V c + dV c at red- 
shift z. The deceleration parameter q Q is taken to be 0.5 
. We have chosen flat models for the universe ( q o =0.5) 
which are popular (Ostriker & Steinhardt, 1995) as they 
can be produced naturally in any inflationary scenario and 
require no fine tuning of parameters. These models are 
also consistent with independent constraints from various 
cosmological observations, e is the fraction of halos that 
give rise to absorption lines; for the galactic halos it is the 
fraction of halos having sufficient gas to produce absorp- 
tion lines. For these, e is taken to be 0.69, which is the 
fraction of late type galaxies (Postman & Geller 1984). 
For mini-halos e = 1 . The redshift distribution of lines of 
a particular ion with column density greater than Ni on is 
given by 



dN(z,N ion ) 
dz 



H 



(1 + z) 



1/2 h 

v. 



(V C) z) x 



1 

J dp A(N ion ,M) edV c . 



(14) 



We assume two models for the velocity distribution of ha- 
los : 

(i) Mass condensations in the CDM universe obtained 
using the Press-Schechter theory (hereafter PS distribu- 
tion) : From the CDM models of structure formation (Mo, 
Miralda-Escude & Rees 1993) the mass and circular ve- 
locity of a halo are related to the comoving radius r Q and 
redshift z by, 

M = p r Q , V c = 1.67 (1 + z) 1 / 2 Ho r Q , (15) 



where p is the mean density of the universe at the present 
time. The number of halos per unit volume with circular 
velocity between V c and V c + dV c is given by 



n(V c ,z)dV c 



-3(1.67) 3 £ c H Q (l + z) 5 / 2 din A 
(27r)3/ 2 V4A(r ) d In V c 

-S 2 {\ 4- y) 2 
x cxp ( ^L±_L) dV c . (16) 

Here S c = 1.68 and the functional form of A(r Q ) for the 
CDM power spectrum of density perturbation is 
A(r Q ) = 16.3 bg 1 (1 - 0.3909 r" 1 + 0.4814 r^' 2 )" 10 , (17) 
where b g is the bias parameter which is taken to be 1 in 
our calculations. 

(ii) Observed (Schechter) luminosity function of galaxies 
(hereafter LF distribution) : The number of spiral galaxies 
per unit volume, with luminosity between L and L+dL is 
given by, 

n(L)dL = </>*(^) a exp(--^)dL. (18) 

Here <jf and a give the normalization and the index of 
faint end slope of the luminosity distribution. L* is the lu- 
minosity of the galaxies having circular velocity V* which 
we take as 250 km s _1 . The luminosity L is related to 
circular velocity V c as = (^-) 2,6 for spiral galaxies. 
Bromley et al. (1998) have determined the values of the 
luminosity function parameters for six spectral type of 
galaxies in two different environments. These values are 
used here to calculate the total number of galaxies as a 
function of the circular velocity and redshift. 

It is possible that the number density of galaxies has 
been higher in the past and reduced with time due to 
merger of galaxies. We also consider this possibility. The 
number density of galaxies after merging is given by n(St) 
= f(<5t)n Q (Park & Gott 1997) where n Q is the present 
number density of galaxies and St is the look back time of 
the universe given by, St = |(1 — (1 + z) -1 - 5 ). We assume 
that there is no evolution of the velocity dispersion. The 
strength and time dependence of merging is described by 
the function f(<5t) which is assumed to be exp(Q(<5t)) . 
The merging rate Q roughly corresponds to the average 
number of pieces at z ^0.65 that merged to produce one 
galaxy at z = 0. The value of Q varies between 2 to 4 for 
different models of the universe; here we take Q=3.5 as 
that gives the best match between model predictions and 
observations. 



4. Constraining the model parameters 

The free parameters in our model are rj c \(0), the number 
of clouds per unit volume at the centre, thickness of the 
disc, t and a, which determines the rate of evolution of 
the disc size. Increasing t decreases the particle density 
and therefore decreases the fraction of neutral hydrogen. 
A given value of Nhi is thus produced at smaller values 
of the impact parameters for larger values of t reducing 
the number of DLAs per unit redshift. Change in a on the 
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other hand changes the slope of the redshift distribution 
of the DLAs. 

The observed number density of DLA absorbers as a 
function of redshift has been obtained by Lanzetta et al. 
(1995) and Storrie-Lombardi et al. (1996). The number 
density evolution of DLASs with Nhi > 20.3 for our model 
is shown in Fig.l for 4 values of a and three values of t 
for PS distribution and a=0.9 and t=0.01 Rd for LF dis- 
tribution. We find that the model results for the redshift 
distribution of DLAs match best with the observations for 
rcdshifts smaller than 4 for t ~ 0.01 Rd and a ~ 0.9. These 
are the values we use in the rest of the calculations. We 
however, note that the observed number of DLASs at z 
~ 4 is much larger than the prediction of all our models. 
Note that as the DLASs are formed purely by the disc, 
their number does not depend on ?7 c i(0). 

Storrie - Lombardi et al. (1994) have studied that 
the redshift distribution of LLSs and have found it to 
be roughly a power law with slope 1.55 over the redshift 
range 0.40 to 4.69. However, they find that the distribu- 
tion becomes steeper for redshifts z > 2.5. Their values 
are plotted in Fig. 2. The figure also shows the results of 
our calculations for the number density evolution of LLSs 
for two values of f (which determines the number of clouds 
per unit volume at the centre of the galaxy, see eqn 2.7) 
and for the halo radius Rh = 1000 and 500 kpc for PS dis- 
tribution. The radius of the halo does not make any differ- 
ence to the distribution. This is because the contribution 
from the outer halo to HI column density is negligible for 
los producing LLSs. The value of f=0.0035 roughly agrees 
with the observed values for z<4 for PS distribution. The 
values for LF distribution are also plotted in the figure 
for f=0.0015 and Rh=500 kpc which are in reasonable 
agreement with the observed values. We use f= 0.0035 for 
PS distribution and f=0.0015 for LF distribution in what 
follows. Here again we note that for z^ 4 the observed 
number of LLS is considerably larger than the model pre- 
dictions. For f = 0.0035 we estimate that a los with impact 
parameter < 100 kpc through halos with V c > 150 km s _1 
will have > 90% probability of intersecting a cloud. The 
probability is however less than 90 % for f =0.0015. The 
value of f obtained for PS distribution is therefore consis- 
tent with the lensing data. 

The number density evolution of LAASs is displayed 
in Fig 3 for Rh = 1000 and 500 kpc. As seen, this is sen- 
sitive to the value of halo radius as los producing LAASs 
does get significant contribution from the outer parts of 
the halo. The figure also shows the observed redshift dis- 
tribution (taken from Kim et al. 1997) for column density 
> 10 1377 cm~ 2 , which is also the lower limit used for the 
theoretical curves. The distribution is consistent with the 
observed distribution for Rh =500 kpc, for both PS and 
LF distribution. We, however, note that the contribution 
of the minihalos is independent of the value of their radius 
as long as it is larger than 250 kpc. So we have assumed 



the minihalo radius to be 250 kpc and the galactic halo 
radius to be 500 kpc in the rest of the calculations. 

It can be seen from the figures 1, 2 and 3 that the 
use of LF distribution without merger is inconsistent with 
the observations of the redshift evolution of DLASs, LLSs 
and LAASs. The observations thus do indicate a strong 
evolution in the comoving number density of galaxies. 

The column density distribution function for hydrogen, 
f(Nni), is defined as the number of absorbing systems per 
unit column density per unit redshift path which is de- 
fined by X(z) = §[(1 + z) 3 / 2 - 1] for q = 0.5. A com- 
parison of the observed column density distribution with 
the predictions of our model is presented in Fig. 4 for the 
column density range 10 12 — 10 22 cm~ 2 at redshift z = 
2.5. The observed data are from Storrie-Lombardi et al, 
(2000). A reasonably good match between the theoretical 
and observed distributions is obtained for both PS and LF 
distribution. 

5. Redshift distribution of heavy element 
absorption systems 

For determining the redshift distribution of HEASs, one 
has to consider the change in abundance of heavy ele- 
ments with redshift. It seems possible that the abundance 
in galactic halos has been increasing continuously with 
time because of in situ star formation (Khare & Rana 
1993; Chiba & Nath 1997). We therefore assume the abun- 
dance of heavy elements to depend on redshift as Z(z) = 
Z(0)(l+z) _(5 in the galactic halos as well as galactic discs 
and consider several values for 5. We have also consid- 
ered another case of an assumed abundance gradient Z(r, 
z) = Z(0, z)e~ 2r / Vc for the heavy elements produced in 
situ in the galaxies. As there is no star formation in the 
minihalos, the abundance in these halos has been assumed 
to be [Z/H]=-2.5 which has been observed in low column 
density Ly alpha clouds (Songaila & Cowie, 1996). 

The redshift evolution of CIV for N C i V > 10 13 c m - 2 
is shown in Fig. 5 for 5 = and 3 and for 5 = 3, 4 with 
an assumed abundance gradient for PS distribution. For 
LF distribution lines are shown for S — and 3 with an 
abundance gradient without merger and 5 — 3 and 4 with 
an abundance gradient with merger. In order to obtain 
the value of dN/dz from observations, a maximum likeli- 
hood analysis was performed with the data collected from 
the literature (references given in Das & Khare 1999). 
The data include components of DLASs also. These are 
included in the data, since most of the components are 
likely to arise in the galactic halos. For this analysis we 
have only used the absorption systems having a relative 
velocity w.r.t. the QSOs larger than 5000 km s -1 . We have 
thus excluded any intrinsic absorption systems that may 
be present in the QSO spectra. The same holds for any 
other observational data used in this paper. The observed 
values are plotted in the Fig. 5. Songaila (1998) has re- 
ported observations of CIV toward 13 QSOs. The value 
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of dN/dz calculated for her data by maximum likelihood 
analysis is also plotted in the figure. As seen from the fig- 
ure, high S values (3,4) are indicated by the observations. 

We have constructed a sample of observed Mgll col- 
umn densities in individual components in the absorption 
systems at redshift z^0.5 from the literature. We have 
thus included all the observations of Lanzetta & Bowen 
(1992); Carswell et al. (1991) and Pctitjean & Bergeron 
(1990) and performed a maximum likelihood analysis. The 
resulting number of systems per unit redshift for lines with 
Njvigii > 10 1315 cm -2 seems to require S > 0. However, the 
observed value is for z ~ 1. More data at higher redshifts 
are required to draw any firm conclusion about the value 
of S. 

6. Column Density Distribution of heavy element 
lines 

The density distribution function for Mgll, fN Mg m i s 
shown in Fig. 6, for the two velocity distribution of the 
absorbers. Observed values which are the result of maxi- 
mum likelihood analysis of the data sample collected from 
Lanzetta & Bowen (1992), Carswell et al. (1991) and Pe- 
titjean & Bergeron (1990) are also plotted in the figure. As 
evident from the figure, there is a reasonably good agree- 
ment of the model results with the observations with 6=3. 
However the column density distribution, being plotted on 
logarithmic scale is not very sensitive to the delta values 
used. 

The column density distribution of CIV lines for 2.52< 
z < 3.78 has recently been obtained from high resolution 
Keck observations by Songaila (1998). This distribution is 
roughly a power law with a slope of -1.5 for column den- 
sities larger than 6 x 10 12 cm~ 2 for which her sample is 
complete. We have plotted model results for the column 
density distribution at z=3 for several 5 values in Fig. 7 for 
both PS and LF distribution. Here, as in case of redshift 
evolution, results for high S values i.e. 3 and 4 are in better 
agreement with the observations. However, at lower col- 
umn densities the model predicts many more lines than 
the observed number. The observed data, however, may 
be incomplete below 6 x 10 12 cm~ 2 as noted by Songaila 
(1998) and it is possible that the number of small column 
density lines may actually be considerably larger. An in- 
crease in the rate of abundance evolution does not serve 
the purpose, since it reduces the number of higher column 
density lines more than that of lower column density lines. 

7. Conclusion 

Pure halo or pure disc will not produce sufficient variety 
in the properties of the QSO absorption lines to explain 
the ensemble of observed lines. To better understand the 
structure of the absorbers and to set constraints on the 
their properties we have presented a model which includes 
minihalo, galactic halos and galactic discs. We compared 



various observed distributions of the absorption lines with 
the predictions from our model. 

The main result of the work presented in this paper is 
that the statistical distributions of various classes of ab- 
sorption systems can be reasonably well explained by the 
three component models which are expected to arise natu- 
rally in CDM cosmology. We do not claim the model to be 
unique and there could very well be other sets of parame- 
ters which can also explain the observations equally well. 
We also note that some of the statistical properties like the 
redshift distribution of Mgll lines are not very well deter- 
mined and may change when more observations become 
available. However at present we find the observations to 
be consistent with the CDM cosmology. 

Some of more quantitative results are as follows. 

(1) The ratio of the disc radius to virial radius decreases 
with time roughly as (1+z) 9 . 

(2) For our assumed baryonic mass of the halo clouds 
(<~ 10 6 M Q ) there are roughly up to a few hundred clouds 
in the halos. Thus roughly a few percent of the baryonic 
halo mass is in the form of a clouds. 

(3) The Lyman alpha lines with column density 
> 10 13 - 77 cm~ 2 seem to require the radius of the galactic 
halos to be around 500 kpc assuming this radius to be 
independent of redshift. 

(4) The column density distribution of HI and Mgll 
is rather insensitive to variation in various parameters 
and is in an excellent agreement with most models, The 
agreement of the observed column density distribution of 
HI with the model results for over 8 orders of magnitude 
in column density does strengthen the faith in a unified 
model like the one presented in this paper. 

(5) The redshift distribution and column density distribu- 
tion of CIV lines does require strong metallicity evolution 
with redshift as (l+z)~ 5 , S being > 3, along with an 
assumption of the presence of an abundance gradient in 
the halo, the abundance decreasing exponentially with 
distance. However, we note that CIV lines mostly arise 
in clouds embedded in the halos and give us information 
regarding the abundance in the halo and not in the discs. 
The observations of Mgll lines which get contribution 
from the disc also are not sufficiently well known to 
give direct evidence for or against such an abundance 
evolution in the discs. 

(6) Observed luminosity function of galaxies is consistent 
with the observations of QSO absorption lines only if one 
assumes that the comoving number density of galaxies 
was larger in the past and reduced with time due to 
merging. 



8. Discussion 

We have presented above a unified model for QSO ab- 
sorbers where various classes of QSO absorption systems 
are produced (apart from the few low column density Ly 
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a lines produced by the minihalos) by lines of sight cross- 
ing spiral galaxies at different values of impact parame- 
ters. As seen above the model is able to reproduce the ob- 
served redshift and column density distributions of various 
species of the absorption systems reasonably well. Several 
doubts have been expressed over such a unified scheme 
in the literature. Bowen et al. (f 995) noticed considerable 
effect of environment on the absorption lines produced 
by galaxies. They, therefore cautioned against construct- 
ing models of generic galactic halos. Diverse morphologies 
have also been observed for DLA galaxies (Le Brun et al. 
1997; Turnshek et al. 2000; Kulkarni et al. 2000). 

LLS are known to originate in bright galaxies (Stei- 
del et al. 1995). However, it has been suggested lately 
that DLASs are produced by low surface brightness galax- 
ies. High redshift DLASs have not been detected in Ly 
a or H a emission except in a few cases, while low red- 
shift DLASs have been observed to be associated with low 
surface brightness galaxies (Turnshek et al. 2000, Pettini 
et al. 2000). Rao and Turnshek (2000) found no evolu- 
tion in the HI content of DLASs over the redshift range 
of 4.0 to 1.0. As a lot of star formation has been pre- 
sumably taking place in the bright spiral galaxies over 
this period it would have led to a significant decrease 
in HI content of these galaxies. They therefore suggest 
that DLASs and star forming galaxies may be different 
populations. Pettini et al. (1999) showed that DLASs do 
not show significant chemical evolution from z~3 to 0.3. 
They also therefore suggest the DLASs to be different 
from bright galaxies which undergo extensive star forma- 
tion and evolution. They suggest that the bright galaxies 
are under represented in DLA samples precisely for this 
reason. The bright galaxies having higher metallicity and 
therefore higher dust content will render the QSOs ex- 
tinct. Savaglio (2000) noted that the metallicity discussed 
by Pettini et al. (1999) (being weighted by column den- 
sity) is dominated by high column density systems which 
may not be uniformly distributed over the entire redshift 
range. Restricting to a sample of 50 low column density 
(N H i < 10 20 8 cm~ 2 ) DLASs which are presumably dust 
free she found evidence for chemical evolution. The total 
sample of 75 QSOs also shows evolution if dust depletion 
is corrected for. She found the global metallicity to change 
from Z Q /30 to Z Q from z ~ 4.1 to 0.4. It is therefore possi- 
ble that the observed DLA sample is biased towards LSB 
galaxies. 

Theoretical studies (Mo, Mao & White 1998; Mo & 
Mao 1999) show that the halo spin parameter plays an 
important role in determining the properties of galaxies 
forming within halos of dark matter. Halos of small angu- 
lar momentum give rise to compact high density systems 
while halos of low mass and high angular momentum form 
discs with a low surface density of baryons. The former 
are sites of star formation while the later dominate the 
absorption cross section . 



Our model does not take into account the presence 
of dust and assumes the chemical abundance to be inde- 
pendent of the mass of the galaxies at a given redshift, 
though in some of the models we do assume the chemical 
abundance to decrease with radial distance, the rate of 
decrease being higher for halos with lower circular veloci- 
ties. However, even if high metallicity systems are indeed 
under represented in the observed DLA sample the rate 
of chemical evolution deduced from our model will not 
change. This is because the rate of chemical evolution ob- 
tained here is from the redshift evolution of the CIV lines 
with column density > 10 13 cm -2 which is decided by 
the los passing only through the halo. Thus the chemical 
evolution deduced in our models is confined to the galac- 
tic halos and does not necessarily reflect a corresponding 
chemical evolution in the discs. The thickness of the discs 
may, however, have to be decreased further in order to 
account for the DLASs which arc unobserved due to dust 
extinction. 
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Figure caption: 

Fig.l: Redshift distribution of DLASs with column 
density greater than 10 20 3 cm~ 2 . Dashed lines marked 
1,2,3,4 arc for a=0.0,0.7,0.8,0.9 respectively for PS 
distribution. Upper and lower dotted lines are for t = 
0.005 R<2 and 0.02 R^ respectively for the PS distribution. 
All other lines are for t = 0.01 R^. Upper and lower solid 
lines are for LF distribution with and without mergers 
respectively. Points marked 'x' are observed points from 
Storrie - Lombardi et al (1996). 

Fig. 2: Redshift distribution of LLSs. Lower and upper 
solid lines are for the PS distribution for Rh=500 kpc 
for f=0.0030 and 0.0035 respectively. Upper and lower 
dash-dotted lines are for LF distribution and Rh=500 kpc 
with and without mergers respectively with f=0.0015. 
Points marked 'x' are observed points from Storrie - 
Lombardi et al (1994). 

Fig. 3: Redshift distribution of LAASs with column 
density greater than 10 13 77 cm -2 . The solid line is the 
observed distribution taken from Kim et al (1998). Upper 
and lower dotted lines are for Rh=1000 and 500 kpc 
respectively with f=0.0035 for PS distribution. Upper and 



lower dashed lines are for Rh=500 for LF distribution 
with and without mergers respectively with f=0.0015 
Fig. 4: Column density distribution of HI lines at z=2.5. 
Points marked 'x' are observed data taken from Storrie- 
Lombardi et al. (2000) and the dotted line is the best fit 
to the points. The solid line is for PS distribution with 
a=0.9, f=0.0035, Rh=500 kpc. Upper and lower dashed 
lines are for LF distribution with and without mergers 
respectively with a=0.9, f=0.0015, R H =500 kpc. 
Fig. 5: Redshift distribution for C IV lines with 
Nciv > 10 13 cm~ 2 . Upper and lower long dashed lines 
are for 6=0 and 3 respectively for PS distribution. Upper 
and lower short dashed lines are for 6 =3 and 4 with 
an abundance gradient respectively for PS distribution. 
Upper and lower dotted lines are for 6 = 3 and 4 with 
an abundance gradient respectively for LF distribution 
with mergers. Upper and lower solid lines are for 5 = 
and 3 with an abundance gradient respectively for LF 
distribution without mergers. The triangle shows the ob- 
served value for the data set collected from the literature 
as described in the text, and the square represents the 
value for the data from Songaila (1998). 
Fig. 6: Column density distribution for Mgll lines at z 
= 0.5. Solid line is the observed distribution which is 
the result of the maximum likelihood analysis of the 
data taken from the literature. Upper and lower dashed 
lines arc for 6= 3 and 6= 3 with abundance gradient 
respectively for PS distribution. Upper and lower dotted 
lines are for 5 = and 3 respectively for LF distribution 
without mergers. Long dashed and dash dotted lines are 
for 5 = and 3 respectively for LF distribution with 
mergers. 

Fig. 7: Column density distribution for C IV lines at z 
= 3.0. Long-dashed line is for 6 = 4 and short-dashed 
line is for 6=3 with an abundance gradient for PS 
distribution. Dotted and long-short dashed line are for 
5 = 3 and 4 respectively with an abundance gradient 
for LF distribution with mergers. Solid line is for 6 = 
for LF distribution without mergers. Short dash-dotted 
and long dash dotted lines are for 5 = 3, 4 respectively 
with an abundance gradient for LF distribution without 
mergers. Triangles and squares are the observed values 
from Songaila (1998) for z< 3 and z > 3 respectively. 
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